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Abstract 

The properties of low-frequency convective fluctuations and transport are investigated 
for the boundary region of magnetized plasmas. We employ a two-dimensional fluid model 
for the evolution of the global plasma quantities in a geometry and with parameters relevant 
to the scrape-off layer of confined toroidal plasmas. Strongly intermittent plasma transport 
is regulated by self-consistently generated sheared poloidal flows and is mediated by bursty 
ejection of particles and heat from the bulk plasma in the form of blobs. Coarse grained 
probe signals reveal a highly skewed and flat distribution on short time scales, but tends 
towards a normal distribution at large time scales. Conditionally averaged signals are in 
perfect agreement with experimental measurements. 

It is well established that the cross field transport of particles and heat near the edge of 
magnetically confined plasmas is strongly intermittent. This is observed in a variety of devices 
including linear HI 13 as well as toroidal configurations EH). Detailed investigations of the spa- 
tial fluctuation structure have revealed strong indications that the intermittent nature of particle 
and heat transport is caused by localized structures in the form of plasma "blobs" propagating 
radially far into the scrape-off-layer (SOL) of toroidal plasmas [ 2', '5', '6'! . It was suggested that 
this is caused by a dipolar vorticity field formed by the charge separation in a density blob due 
to guiding-center drifts in a curved inhomogeneous magnetic field f?]. 

In this contribution we will provide a self-consistent description of the intermittent particle 
and heat flux and link it to the of the emergence and evolution of such blob-like structures. We 
base our investigations on a novel model for interchange turbulence in slab geometry for the 
outboard midplane of a toroidal device [8|. The model includes the self-consistent evolution 
of the full profiles in the edge/SOL. A local Boussinesq-like model, where the "background" 
profile is separated from the fluctuations, fails to provide a realistic description. The geometry 
comprises distinct production and loss regions, corresponding to the edge and SOL of magne- 
tized plasmas. The separation of these two regions defines an effective last closed flux surface 
(LCFS), though we do not include magnetic shear in our model. In the edge region, strong pres- 
sure gradients maintain a state of turbulent convection. A self-regulation mechanism involving 
differential rotation leads to a repetitive expulsion of hot plasma into the SOL, resulting in a 
strongly intermittent transport of density and heat. 
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Figure 1 : Geometry of the simulation domain showing the forcing region to the left, correspond- 
ing to the edge plasma, and the parallel loss region to the right, corresponding to the scrape-off 
layer. Parameters are 5 = 8 and ^ = 1 . Data time series are collected at the probe positions Pi. 



The model derives from the continuity equations for the electrons, the electron temperature 
and the quasi-neutrality condition. Assuming cold ions and neglecting electron inertia effects, 
we obtain a three-field model for electrostatic perturbations of the full particle density n, electric 
potential (|) and electron temperature T . Using slab coordinates with z along the magnetic field, 
X in the radial and y in the poloidal direction we obtain |i8J, 
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where time is normalized by l/cOc/ and spatial scales by = Cs/cDd- The density n and tem- 
perature T are normalized to fixed characteristic values at the outer wall. We further define the 
advective derivative, the magnetic field curvature operator and the toroidal magnetic field by 
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The vorticity Q. = V^(|), the inverse aspect ratio e = g/Rq and C, = Ps/Ro where a and Rq are the 
minor and major radius of the device. The terms on the right hand side of the equations describe 
external sources S, parallel losses along open field lines through the damping rates o lIHl, and 
coUisional diffusion with coefficients v. The geometry and boundary conditions are sketched in 
Fig.m 

In the absence of external forcing and dissipative processes the model equations non-linearly 
conserves the global energy to lowest order in C, 
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where the integral extends over the whole plasma layer. We define the kinetic energy of the 
fluctuating and mean components of the flows, 



where the zero index denotes an average over the periodic direction y and the spatial fluctuation 
about this mean is indicated by a tilde. The linearly damped mean flows, vq = d^o/dx, does not 
yield any radial convective transport and hence form a benign path for fluctuation energy. The 
energy transfer rates from thermal energy to the fluctuating motions, and from the fluctuating to 
the mean flows, are given respectively by 



Note that Fp is directly proportional to the domain integrated convective thermal energy transport, 
while Fy shows that structures tilted such as to transport positive poloidal momentum up the 
gradient of a sheared flow will sustain the flow against coUisional dissipation ll9l fT0lfTT1l . 

In the following we present results from numerical simulations of the interchange model 
using parameters relevant for SOL plasmas. = 2Ly = 200 and the LCFS is located at Xlcfs = 
50. The parameters are £ = 0.25, ^ = 5 x 10~^, and v = 10^^ for all fields. The parallel loss 
rate of temperature is assumed to be five times larger than that on the density and vorticity, 

= <^Q. = = 3C,/2nq, since primarily hot electrons are lost through the end sheaths. a„ 
and Of! correspond to losses at the nearest target plate over a distance of L|| = 2KRoq/3 (one third 
of the connection length) with the acoustic speed Cs, where ^ = 3 is the safety factor at the edge. 
Finally, the radial line-integral of the sources S„ and St equals 0.1. For the numerical solution 
the spatial resolution is 512 x 256 grid points in the radial and poloidal directions, and the time 
span of the simulation is 2 x 10^. 

We have performed several runs with varying parameters, showing that the qualitative be- 
havior is robust, whereas the quantitative results depend on the parameters and particular on the 
value of the coUisional diffusivities. The general observation is that the turbulent flux is strongly 
intermittent: quite periods are interrupted by strong bursts of particle and heat fluxes. This is 
correlated with the kinetics in the fluctuations, as is shown in Fig. El We observe that the convec- 
tive energy and thermal transport appears as bursts during which particles and heat are lost from 
the edge into the SOL region. As discussed in Refs. ||9l[T0l[TTl, this global dynamics is caused 
by a self -regulation mechanism in which kinetic energy is transfered from the fluctuating to the 
mean components of the flows, and subsequently damped by coUisional dissipation. The thermal 
energy ejected in a bursty manner from the edge and into the SOL region, will eventually be lost 
by transport along open field lines. The characteristics time between the bursts is related to the 
viscous diffusion (compare with Fig. 3 in Ref. |8|, where the value of v is 5 x 10^^). We further 
verified that the self-sustained poloidal flow profiles are strongly sheared in the edge region, and 
have larger amplitudes during the strong fluctuation period. 

The statistics of single-point recordings at different radial positions Pi indicated in Fig. [T] 
agree very well with experimental measurements. In Fig.|3we present the probability distribution 
functions (histogram of counts) (PDF) of the density signals taken from a long-run simulation 
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Figure 2: Evolution of the kinetic energy contained by the mean U and fluctuating K motions 
[Eq. ([T])] and the collective energy transfer terms [Eq. Q] . 

containing more than a hundred strong burst events. It is notably that the PDF at the first probe 
inside the LCFS is close to a Gaussian with skewness 0. 12 and flatness factor 2.97, while the PDF 
becomes more flat and skewed further into the SOL. This indicates the high probability of large 
positive fluctuations corresponding to blobs of excess plasma. The skewness and flatness factors 
grow through out the SOL and take values up to 4 and 25, respectively. The PDF's in the SOL 
have similar structure with a pronounced exponential tail towards large values, a characteristic 
feature of turbulent convection in the presence of sheared flows [9 , 10|. 

We have also considered the coarse-grained PDF which is obtained by averaging the signal 
over time intervals of lengths x and constructing new time records with a time resolution of x: 
riT^it) = (1/x) //_X/2^"(^') ■ coarse grained PDF's (PDE^) for the signal at P^, are also plot- 
ted in Fig.|3lfor increasing values of x. We observe that PDFx approaches a Gaussian distribution 
when X is exceeding the averaged time interval between bursts, which is roughly 10^. This shows 
the absence of self- similarity, which is characteristic for an intermittent signal (see, e.g., 1 12 1). 

The conditionally averaged temporal wave forms of the density calculated from the same 
signals and the radial velocity field v^, using the trigger condition n — n> An^^^ at each individual 
point, are presented in Fig.|4| For the density signal an asymmetric wave form with a sharp rise 
and a relatively slow decay is clearly seen, as also observed in experimental measurements O 
|510|. The maximum density excursions significantly exceed the background level, and decay 
as the structures propagate through the SOL. By using a negative amplitude for the conditional 
averaging very few realizations results, confirming the presence of blob-like structures. For 
the velocity signal we observe that the radial velocity is positive (directed radially outwards) 
in the blob. In the edge region it takes weak negative values both before and after the blob. 
Also this result agree with experimental observations |5 1. We note that the maximum value of 
decreases on passing the LCFS and then increases to a maximum value of 0.046 at Pt,, after which 
it slowly decays. From two-dimensional animations we clearly observe the radial propagation of 
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Figure 3: In the left panel is shown the probability distribution functions of particle density 
measured at seven different radial positions Pi as shown in Fig. [T] To the right is the coarse 
grained probability distribution function taken at probe 3, averaged over the time intervals x 
indicated by the labels. With increasing x the skewness decreases as: 2.6, 1.9, 0.73, 0.16, and 
the flatness factor decreases as: 12.0, 7.8, 3.9, 3.2. For both plots the vertical axis shows count 
numbers on a logaritmic scale, n designates the averaged density. 



blob-like structures for the density and temperature fields, while the vorticity displays a dipolar 
structure as expected from theory and experiment (cf. Ref. |8|). From such animations and 
radial correlations we find that the radial propagation velocity of the blob structures corresponds 
to around O.OSc^ consistent with Fig. HI but with a large statistical variance in agreement with 
experimental measurements 013151. 

By combining the conditional evolution of n and Vx in Fig. |4l we deduce that the blobs are 
carrying a large particle flux. We have examined the PDF of the particle flux averaged over the 
periodic y-direction (the flux surface) at different radial positions. The PDF's are quite similar 
and strongly skewed with a flat exponential tail towards positive flux events, showing that the 
flux is dominated by strong bursts. The tail of the PDF was found to be well fitted by an extreme 
value distribution II13I . By coarse graining the PDF as described above we observe a similar 
behavior as for the local density fluctuations: the distribution approaches a Gaussian for large 
time scales. 

We have demonstrated that a two-dimensional model for interchange turbulence provide re- 
sults in good agreement with that reported from experimental investigations of SOL turbulence 
and transient transport events [T TrSI. An important feature of the model is the spatial separation 
between forcing and damping regions. Our results are in quantitative agreement with experimen- 
tal measurements of field-aligned blob-like structures propagating far into the scrape-off layer. 
The associated intermittent transport may have severe consequences for magnetic confinement 
experiments by producing large heat bursts on plasma facing components. 
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Figure 4: Conditionally averaged wave forms of the particle density (left panel) and the radial 
velocity measured at seven different radial positions Pi as shown in Fig.[Tl using the condition 

n{xp.) -n{xp) > 4ntms{xpi). 
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